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SUMMARY 


In previous reports By the author, the method of operators has been 
applied to the investigation of two-dimensional irrotational flow patterns 
of an ideal compressible fluid, mostly in the subsonic case. Using a 
similar approach in the present paper, the author defines two different 
operators which generate supersonic two-dimensional flow patterns from 
differentiable functions of one real variable. (An operator is any rule 
by means of which one function is converted into another. This concept 
is discussed in some detail in the body of this paper.) Since operators 
of the type considered preserve many properties of the functions to which 
they are applied and since the theory of functions of one real variable 
is more extensively developed than that of solutions of the compressi- 
bility equations, the results obtained can be used as a basis for the 
investigation of supersonic flow patterns. 


INTRODUCTION 


By following the line of approach developed in previous publications 
for the study of two-dimensional subsonic flows (see references 1 to 4), 
the method of operators is applied in the present paper for generating 
two-dimensional steady supersonic flow patterns. Using the holograph 
method introduced by Chaplygin, the author considers the stream function 
in the A, 0 — plane, where 


A = h ■'•arc tan [h(M^ — l)-'-/^J — arc tan [(#-l)V2] 
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(see equations ( 5 ) and (6)), M being the local Mach number; 0, the 
angle between the velocity vector and the x-axis; and h, a constant 
characteristic of the fluid. 1 

The stream functions, considered as a function of A and 0, 
satisfy a linear differential equation 


*aa - *00 + 4n i*a = 0 


Two operators which generate solutions of this equation are defined 
in the present paper. The first operator is obtained by using Eiemann's 
function. (See section entitled "Operator Obtained by the Use of 
Riemann’s Function.") The second operator is derived by means of the 
theory of integral operators. In the section dealing with the second 

operator, a set of functions H(A) , E m I1 (A) , n = 1, 2, . . . , are 
determined such that 



<*• 

00 


P(f) = lim «< 

H(A) 


> 

m — » 00 


n=l 



1 At every point (x,y) of a steady two-dimensional flow the 
velocity vector, the Cartesian components of which are q cos 0 
and q sin 0, is defined. If \| r(x,y) denotes the stream function 
8111 p( x *y)> the density, the coordinates x and y can be expressed 
as functions of q and 9 (and hence, using relations (5) and (6) 
of A and 0) by solving the equations 


and 


q cos 0 


-1 cty 
dy 


q sin 0 



Substituting these expressions for x and y into i(x,y), the 
function ^(A,0) = i [x(a, 0), y(A,0)] is obtained, i.e., the stream function 
in the A, 0— plane. 
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where 


f M (t) 


n 

Ja 


f^(S) = f(l) 


I = A + 0 or | = A — 0 


(where f is an arbitrary differentiable function of one real variable £ ) 
represents a solution of the equation 


*AA “ *00 + ^ N 1*A = 0 


so that P(f) can be interpreted as the stream function of a supersonic 
flow. The equation of a streamline of the flow in the physical plane is 
then given in parametric form: 


x 



dA 1 . 

cos 0 V, 

cLq pq 


0 


sin 0 


) 


d0 


(m 2 - l) 


L pi 


cos 0 ^ — sin 0 

e dA pq A 


dA 


V \ 

y = ( - t* 3 in 0 — + — cos 0 ) d0 
\P A dq pq *0 / 


(M 2 - l) ' ■ dq 1 , 

sin 0 — + — cos 0 


pq 2 ' w . cLA.: Pq " A 


dA 


where \|r(A,0) = c = Constant. The representation P(f) is investigated 
in the present paper, and, in particular, methods for .the actual 
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determination of E^ 11 ^ are discussed as well as the domain in the 
A, 6— plane where the order of summation and passage to the limit can he 
inverted . 

The evaluation of the operators , that is, the determination of the 
flow pattern when the function f is given, involves lengthy computa- 
tions. However, in a manner similar to that described in reference 4 
for the subsonic case, it is possible to prepare, once and for all, 
tables which greatly facilitate numerical computation of flow patterns. 

The results obtained in the present paper can be used as a basis 
for the investigation of supersonic flows. 

The next problem which arises in this connection consists in deter- 
mining for a desired flow the functions of one variable which have to be 
inserted into the operator in order to obtain the desired boundary shape 
or to decide that such a flow does not exist. The author expects to 
treat these problems in subsequent publications. See also the SUMMARY 
REMARKS presented at the end of the present paper. 

This investigation was conducted at Brown University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

The author would like to express his sincere appreciation for the 
assistance and advice he received from Mr. Bernard Epstein and to thank 
Mr. tfeurice H. Slud for his interest and valuable help in connection with 
the present investigation. 


SYMB0I5 


a 


velocity of sound in compressible fluid 


a o 


velocity of sound at a stagnation point 


h = 


\ 


k - 1 
k + 1 


k ratio of specific heat at constant pressure to specific 

heat at constant volume (k = 1.4 for air) 


<1 

u 

v 




speed of a fluid particle l\lu 2 + r 2 ) 
x— component of velocity vector (q cos 0) 
y— component of velocity vector (q sin 0) 
Cartesian coordinates in physical plane 
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B = \|m 2 - 1 for M > 1 
M Mach number (q/a) 

R(£,q;4 0 ,T| 0 ) Riemann's function for equation (15) 


T = \|1 — M 2 for M < 1 
p = T + iB 


Z = X + iA 
(t) 

F x ( 2 a) 
F 2 (2A’ ) 


functions defined by recurrence relations (22) 
see equation (l6) 

see Theorem under section Integral Operator of T!ype Given 
by Formula (ll) 


F 2 ^(2A') 

“n 


£ = a + e 


a sequence of functions converging to F2(2 A')> see equa- 
tion (40) 

coefficients of Taylor series expansion of F 2 (2A' )j see 
equation (39) 


T) = A — 0 


|». = A* + e 

tj' = A* - 6 


\ 


1 

2 



+ i log 
h 



for 


M < 1; log = log e 


A = 7- arc tan (hB) — arc tan B for 
h 


M > 1 



Q 


angle between velocity vector and positive direction of 
x— axis (arc tan v/u) 


potential function 
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\|r stream function 

i|r* "reduced stream function" (see equation (l4)) 

a sequence of functions converging to i|r*; see equation (37) 
density of fluid at a stagnation point 

p density of fluid in motion 

p "reduced density" ^ p / p Q ) 

N^ see equation ( 9 ) 

Note.— Frequently quantities are considered as functions of different 
pairs of variables, so that different symbols should be used to designate 
the functional dependence in the various planes. However, this is not 
done herein, the same symbol being used in each plane. This should cause 
no confusion, for the meaning is clear from the context. 

ANALYSIS 

The Mathematical Problem in Simplified Form in the Case 
of a Supersonic Motion of a Compressible Fluid 

In this report, operator methods are applied to the problem of 
constructing two-dimensional, irrotational, steady flow patterns of a 
compressible fluid at supersonic velocities. 

Just as in the case of an incompressible fluid, the steady, irrota— 
tional flow is completely described by either the velocity potential 0(x,y) 
or the stream function t(x,y) . If it is assumed that the motion is adia- 
batic and that the fluid is a polytropic gas, then the motion is governed 
by the following system of nonlinear partial differential equations: 
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and 


p = p(x,y) 


Here p = p/p Q is the "reduced density" of the fluid, that is, the 
density p divided hy the density p Q at rest; x,y are Cartesian 

coordinates in the physical plane, that is, in the plane where the 
motion takes place; k is the ratio of the specific heat at constant 
pressure to the specific heat at constant volume; and a Q is the 
velocity of sound in the fluid at rest. (See reference 1, equations (24) 
and (25).) 

A fluid motion satisfying equation (l) may show two essentially 
distinct types of "behavior, according as the speed of the fluid is 
greater or less than the speed of sound in the fluid. The mathematical 
reason is readily seen in the following manner. If p and \|r are 
eliminated from equation (l) the result iB a single, nonlinear, partial 
differential equation of the second order for 0: 

(a 2 - + (* 2 - 0/K - 0 < 2 > 


where 


a2 = a o 2 -| (k-l)(0 z 2 + 0 y 2 ) "(3) 

is the square of the local velocity of sound. (See equation (28) of 
reference 1.) The discriminant of equation (2) is 

- ( a2 ~ 0x 2 )(a 2 - 0y 2 ) 

- -a 2 (a 2 - fc 8 - 9 *) - a 2 (q 2 - a 2 ) (4) 

where q^ is the speed at the point x,y under consideration. In regions 
where the velocity is subsonic, that is, q.< a, this equation is of 
elliptic type, whereas in regions where the velocity is supersonic, that 
is, q > a, equation (2) is of hyperbolic type. (See reference 5, pp. 1-4.) 


1 
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The solutions of these two types of differential equation have markedly 
different functional character, and these differences are, of course, 
reflected as the two different types of compressible fluid motion. 

. In spite of these differences, however, the operators which have 
been used for the solution of subsonic flow problems may be adapted to 
the solution of supersonic problems after certain formal modifications. 
Before this is done, however, it is necessary to make a few preliminary 
remarks . 

The differential equations (l), both for subsonic and supersonic 
velocities, were derived under rather restrictive assumptions, namely, 
that the flow is two-dimensional, that the fluid is inviscid, and that 
the motion is irrotational and steady. Bfy these assumptions the problem 
is considerably idealized; but experience has shown that in a great many 
important cases occurring in practice, when the speed is small, the 
consequences of the theory are in satisfactory agreement with observa- 
tions. Moreover, a physical justification for the apparently drastic 
assumption that the motion is laminar is contained in the Yon Mises 
"hydraulic hypothesis." (See reference 6, p. 84.) 

In many instances where the indicated idealization is not a,dmissible, 
the "idealized" theory can at least be used as a first approximation, 
which can be improved by correction terms. So, for instance, the 
boundary— layer theory describes the motion near the boundary where 
viscosity effects cannot be neglected, whereas in the interior of the 
flow pattern the theory of ideal, inviscid fluid can be successfully 
used without any alteration in the case of a comparatively slow motion. 

Unfortunately, when the speed increases considerably, the situation 
changes completely, and the theory of an incompressible fluid can no 
longer be used as a first approximation. Furthermore, at supersonic 
velocities a completely new phenomenon appears, the influence of which 
is not easily assessed. When a compressible fluid flows around an 
immersed body at supersonic velocities, a very characteristic phenomenon 
is the appearance of a shock line. It appears as a curve (or curves) 
which divides the physical plane into two (or more) unconnected infinite 
regions. In the neighborhood of the shock line, the fluid may no longer 
in general be regarded as nonviscous; there is a discontinuity of 
velocity along the shock line; and in the region of the physical plane 
beyond the shock line (except in a few special cases) the flow is no 
longer irrotational. For all these reasons, the differential equa- 
tions (l) may no longer be supposed to represent the actual physical 
situation. Exact equations can be derived, of course, but they are 
extremely complicated nonlinear partial differential equations for which 
no mathematical theory now exists; and the problem presented (which is 
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neither an initial— value problem nor a boundary— value problem, because 
the shock line is not prescribed) cannot be handled by any known mathe- 
matical means. 



It can be expected that in the same manner as the idealized theory 
of an incompressible fluid represents an excellent "first approximation" 
in the sense Just explained for comparatively slow motions, the idealized 
theory of compressible fluid will represent a similar "first approximation 

for motions with considerable speed^ (if necessary, with certain addi- 
tional alterations). 

In this and previous papers on the subject by the author, therefore, 
no attempt is made to find exact solutions of the actual physical problem. 
Instead it is assumed that the part of the flow around the obstacle 
satisfies the conditions under which equations (l) were derived. The 
method of operators is then used to find and investigate exact solutions 
of this idealized problem. The results may be regarded as first approxi- 
mations to the true facts in the sense Just explained. 

2 

It is clear that for many purposes it will be necessary to develop 
a three-dimensional theory. This theory, even in the case of an incom- 
pressible fluid, is only in the very first stages of development. The 
theory of operators and some other modern mathematical means, as the 
author expectB to show in another paper, will make it possible to deal 
with certain problems in the three-dimensional case. These methods, 
besides yielding results for subsonic motions of a compressible fluid 
will also yield results of certain problems (as yet unsolved) for an 
incompressible fluid. 
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The question may also he regarded from another "viewpoint. It would 
often he possible to find a numerical solution of a specific problem, but 
such a result has little theoretical value. A satisfactory theory must 
be capable of doing much more. It must be able to give at least qualita- 
tive information concerning how the motion of the fluid is affected by 
variations in the form of the immersed body, the specific properties of 
the fluid, or any other significant variables; and preferably it should 
also be capable of providing quantitative predictions. The operator 
method of solving differential equations has the advantage that it may 
be used to investigate the functional properties of the solutions it 
yields. (For details, see references 2 and 3*) 

Remark .— From the purely mathematical point of view, an operator is 
a rule which is applied to a function of an appropriate class (us ually 
denoted as class A) to obtain a function belonging to another class 
(class C). A simple example of an operator is the process of obtaining 
a harmonic function of two real variables (class C) by taking the real 
(or imaginary) part of an analytic function of a complex variable 
(class A). This operator preserves many properties of analytic functions 
and is often used as a tool for investigation of harmonic functions. 

A similar relationship may be seen in the operators employed in the 
theory of compressible fluid. These operators appear in the form (equa- 
tion (ll)) for subsonic flows and in a similar form for supersonic flows. 

In the case of subsonic flows, equation (ll) acts on analytic functions 
of a complex variable (class A) and produces complex solutions of equa- 
tion (10) (class C). In the supersonic case, the operator acts on a 
differentiable function of a real variable (class A') and produces solu- 
tions of equation (9) (class C’). These operators are speci al cases of 
operators introduced by the author which convert analytic functions of a 
complex variable into complex solutions of a given linear differential 
equation of elliptic type and differentiable functions of a real variable 
into solutions of a given equation of hyperbolic type. 

The theory of integral operators gives a deep insight into the property 
of the solutions of the compressibility equation. As indicated in 
references 1, 2, and 3, solutions of equation (10) can be written in the 
form 


J ^ E(M,0,t)f ||[i(M) + i0](l-t 2 )J dt/Jl - t 2 


(See equation (55) of reference 1.) Here A.(M) is a function of the 
Mach number M, which in the subsonic case is real, and in the supersonic 

case purely imaginary. The expression E(M,0,t) = E(k,0,t) is a function 
which, in the "simplified case" (Tricomi case), when considered as a 

function of u = t^(X, + i0)/2k, satisfies the hypergeometric equation. 

In the case of the "exact" compressibility equation, behavior and proper- 
ties of E(A,,0,t) are completely analogous to those in the simplified 
case. It is described in sections 4 and 5 of reference 3. 
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The operator method^ has another advantage which is very important 
from the practical point of view. All too frequently the only thing 
which stands in the way of solving an engineering problem mathematically 
is the difficulty of performing the incidental numerical computations . 
Although the use of modern automatic computing machinery in many instances 
reduces considerably the amount of labor needed, it still remains a 
serious obstacle. Whenever possible, the ultimate numerical evaluation 
Bhould be kept in mind when the mathematical analysis is carried out. It 
is a definite advantage when the computations required by a theory may 
be readily performed. It will be seen that when hydrodynamical problems 
are attacked by the operator method, the formulas which occur contain a 
number of standard functions which are the same in all problems of the 
same type. Consequently these functions may be computed and tabulated 
once and for all. When this has been done, the computations required 
for any specific problem are much shorter and simpler. 


The Method of Operators in the Theory 
of Compressible Fluids 

The idea of the method of operators is now explained briefly. 

It has already been stated that the partial differential equations 
of compressibility are nonlinear and that a mathematical theory of such 
equations does not yet exist. Fortunately, it is possible to linearize 
the equations of motion of a compressible fluid by a change of variables. 
Therefore, hydrodynamical problems are not studied in the physical plane 
where the flow occurs but must be investigated in some auxiliary plane 
where the differential equations are easier to handle. This entails 
considerable distortion of the flow pattern and a consequent loss of 
intuitive appeal. 

The process of linearization is conveniently performed in two stages. 
In the first place it is noted that, in general, in a sufficiently emal l 
portion of the flow region, a particle may be Just as unamb iguously 
specified by its velocity vector as by its position, although this will 

The main mathematical idea of the operator method consists in 
employing certain integral operators which tr ans form functions of one 
complex (or real) variable into "classes" of complex solutions of given 
linear partial differential equations and preserves many properties of 
the class of functions to which the operator is applied. Thus, in the 
case of equations of elliptic type, the application of the operators 
results in the derivation of complex solutions of the given solution's 
equation which in many respects behave like analytic functions of a 
complex variable and help in the investigation of real solutions in a 
manner which bears a close analogy to the role of analytic functions of 
a complex variable for real harmonic functions in two variables. 
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not gener ally be true over the flow as a whole* 4 '. The first change of 
variables is, therefore, from the physical plane with coordinates x,y 
to the so-called velocity or hodograph plane in which the polar coordi- 
nates q,0 are used, q being the speed of the particle, and 0, the 
angle between the direction of motion of the particle and the positive 
x— axis. This transformation obviously introduces distortion (see fol- 
lowing figs.) and it is by no means bi-nmique. 


y 



Motion in the physical 
plane. 



The image in the hodograph 
plane of the flow indicated 
in figure at left. 


since it usually happens that the velocity vectors are equal at different 
points of the flow region. A further distortion is now produced by 
plotting log q and 9 as Cartesi an coordinates (logarithmic plane). 
Finally, the substitution 

A =A(q) = | arc tan [h(M 2 - l) 1 ^ 2 ] — arc tan [(m 2 — 1) 1/2 ] (5) 


(see fig. l) where 



In order to make the correspondence between the flow pattern in 
the physical plane and its image in the hodograph plane 1— to-1, it is 
necessary in this case to locate the hodograph on a Riemann surface with 
the proper number of sheets. 
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has the effect of changing the scale along the log q-axis (pseudo- 
logarithmic plane). This introduces an additional distortion hut 
succeeds in reducing the differential equation for the potential func- 
tion to the form 


0 AA “ 000 ~ ^ n i 0 a = 0 ( 7 ) 


where 


Hi = 


k + 1 


M 1 * 


8 (M=-1) 3/2 


The stream function satisfies the differential equation 


*AA - *00 + ^ N 1*A = 0 


( 8 ) 


( 9 ) 


In the subsonic case, a similar procedure leads to the differential 
equation 


*XX + +00 + ^X = 0 


(10) 


(see equation (^7) of reference l), where X = X(q) is a function of 
the speed alone, the form of which is very similar to that of A(q) in 
equation (5). The operator method of finding subsonic compressible flow 
patterns is now easily described. It is based on making suitable use of 
the following formula for generating solutions of equation (10): 


V(X,0) = Tm P(g) 


P(g) = H(2X) 


g ($) 


+ £ Q (n) (2X) 


n 


* — ^ p^n n i 

n=0 ^ n# 


P^n-l 


> c 




U 0 




where £ = X - i0. (See reference 1, formula^ (66).) In this formula 
g represents an arbitrary analytic function of a complex variable £ 

-'Note that in formula (66) of reference 1 , there is a misprint: the 

factor — ^ * should be replaced by * 


2 11 ni 


2 2n n! 
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which is regular at £ = 0; whereas H and the Q n *s are certain 
functions of X only and are the same for all stream functions. They 
may he computed once and for all and the results used in all future 
hydrodynamical problems . (See reference 4.) 

The foregoing formula is employed in the following manner. If any 
analytic function g of a complex variable (regular at £ = 0) is 
substituted into formula (ll), the result is a stream function of a 
subsonic flow pattern. On the other hand, it is well known that an 
analytic function g(log q — i0) may be regarded as defining a flow 
pattern of an incompressible fluid. Consequently formula (ll) may be 
interpreted as an operator which distorts an incompressible— flow pattern 
into a ’ compressible one. (It may be said that the operator acts as a 
distorting mirror which reflects an incompressible— flow pattern as a 
subsonic compressible— flow pattern.) The operator formula is of such a 
form that a knowledge of the incompressible— flow pattern used leads 
directly to both qualitative and quantitative knowledge of the behavior 
of the stream function of the corresponding compressible flow. 

An obvious objection is that, whereas any incompressible flow will 
lead to some compressible flow when a specific obstacle is prescribed, 
there seems to be no indication of how to find precisely the g to be 
substituted into the operator. A procedure to solve this problem 
approximately can be developed (see reference 7), but this is not con- 
sidered herein. However, it is easy to see why the same analytic 
function g(£) which solves the corresponding incompressible problem 
can be expected to yield a first approximation to the compressible 
problem in the case of low speeds (i.e., Bmall Jfech number); for in such 
a case the solution (formula (ll)) is easily shown to yield a flow, very 
similar to the incompressible flow, which agrees with the physical state 
of affairs. 

t 

The operator method as Just described involves the following 
elements : 

(1) It is necessary to have considerable knowledge about 
incompressible— flow patterns, in particular about their Images in the 
hodograph and logarithmic planes. 

(2) The operator formula has to be derived, the functions H 

and have to be determined, and rules must be developed for 

interpreting properties of incompressible— fluid motions as properties 
of the compressible flows given by the operator. 

(3) After a solution of a problem on compressibility has been found 

in terms of the variables (X,0) (i.e., in the so-called pseudo- 

logarithmic plane), means must be furnished for determining the actual 
flow in the physical plane. 

It is obvious that one reason why this method can be so powerful 
in the study of subsonic compressible flows is that the theory of 
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incompressible flows (i.e., the theory of analytic functions of a complex 
variable) is so well known. 

The supersonic case is attacked by perfectly analogous means. As 
has been indicated by the author , there exist various operators for 
transforming solutions of one equation into solutions of another equa- 
tion; and often for different purposes it is convenient to use different 
operators or different forms of the same operator. In the present 
report two operators which transform solutions of the relatively simple 
hyperbolic equation 


dfjr _ _ o 

8 a 2 Se 2 


( 12 ) 


(i.e., a pair of differentiable functions 


g(S) and h(q) 

>• 

where £ = (A + 6 ) and q = (A -0) 


(13) 


of one variable) into solutions of the compressibility equations are 
derived . 


As soon as the function \|r(A ,0) satisfying equation ( 9 ) is obtained, 
the corresponding flow in the physical plane can be obtained exactly in 
the same manner as in the subsonic case. (See reference 1 , section 14.) 

In contrast to the subsonic case, flows which in the logarithmic 
plane satisfy equation (12), A being some function of the speed, have 

not, to the knowledge of the author of this report, been studied. 
Relatively little is known concerning how to choose g anfl h in order 
to obtain in the physical plane a flow around a prescribed obstacle, 
although there is reason to expect that this investigation should be 
simpler than in the subsonic case. However, this question is not con- 
sidered herein, but is postponed to a future publication. 

The following sections of this report are devoted exclusively to 
the determination of two different operators which transform differ- 
entiable functions of one real Variable into solutions of equation ( 9 ). 

In reference 8 Coburn has indicated the extension of the Chaplyginr- 
Karmanr-Ts ien method to the case of supersonic flows. The stream functions 
of such flows satisfy equation ( 9 ). Coburn indicated one characteristic 
property of these flows. On the other hand, the general properties and 
the form of flows (in the physical plane) which are of interest for the 
purpose of the present approach have not been investigated in Coburn’s 
paper . 
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Operator Obtained by the Use of Riemann's Function 

One method of generating stream functions of supersonic flows from 
expression (13) is by the use of Riemann’s function. 

As indicated in the preceding section, the stream function ty(A ,6) 
satisfies equation (9).' If* instead of A and 0, the variables £ 
«nd ^ defined in expression (13). are introduced, and, instead of \|r, 
the "reduced stream function" 


i|r* = i|rexp 


— 

H+ti 

N-Jt) dT 

1 

a 


(HO 


(where a is an arbitrary constant) is considered, then a formal compu- 
tation shows that \|r* satisfies the equation 1 

** + Fi(l + r\H* = 0 (15) 

I 1 ! . 

where 



d(2A) 


k + 1 9(k + 1) 12k 6k - 14 

6k b 6 34 b 2 


+ (4k + 8) - (3k - 1)B 2 


where 

B 2 = M 2 - 1 


(16) 


(17) 


Remark.— Obviously it is allowed to take for the lower limit a 
of the integral in equation (l4) an arbitrary point, since replacing a 
by a*, a* ^ a, means only that \|r* is multiplied by a constant, 

T Pa* 1 


namely 


exp 


N-Jt) dT. 


The most natural choice is to take for a 


Lya J 

the same value as in the subsonic case. (See formula (ill) of reference 1.) 


As is explained in more detail in appendix C, it is possible to 
extend the variable A to the complex values l = X + iA. Then in the 
subsonic case the variable l assumes values on the negative real 
axis X and on the positive imaginary axis A in the supersonic case. 
The integration can be then taken from t = -00. Since N(0) = oo^ the 
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integration must be performed along a path which avoids 2 = 0, for 
example, along the real axis from — « to X^, X^ < 0; then along a 

circle with radius | X^ | and center at 2 = 0; and finally along a 
segment of the positive imaginary axis from t = iJX^J to t = iA. 

Let R(l , t]; I 0 ,t) 0 ) be the Eiemann’s function of equation (15). 

(See, e.g., pp. 311— 317 of reference 9.) Then, according to a classical 
result (p. 316, formula (j*) of reference 9), 


- +*(5 0J n 0 ) B (£,i:5 0 ,i 0 ) 


+ +*y(?o> Y ) 4Y 


+ 



^ # x( X ^o) B(^VSo^o) dX 


(18) 


where 



and so forth. 


Since for \|r*(l 0 ,T) 0 ) can be chosen an arbitrarily prescribed 
constant A and for i|r*(£ 0 ,Y) and i|r*(X,Ti 0 ) can be chosen two arbi- 
trary continuous functions h(X) and g(Y)- ^provided that h(l D ) = g(^ 0 
the formula 


nn 


T|r*(!,n) = AR(^,T);5 0 ,T| 0 ) + 


h(T)B(| 0 ,Y;| 0 ,i, 0 ) dY 




PI 


+ 


J g(*)K(X,ii o ;$o,Tlo) dX 


(19) 


represents an integral operator which transforms two arbitrary continuous 
functions of a real variable into a solution of equation (15). 



18 


WACA TN No. 1875 


(*0.1) 



(^o> ^o) 


(MJ 


•X 


The rectangle within which the values of are determined by 

prescribing its values along the line segments: Y = t) 0 , i Q — X- i 
and X = l Q , ti 0 ^ Y *) . 


As is well known, the Riemann function of equation (15) can he 
written in the form 


00 

= (~l) n 


( 20 ) 


n=0 


where 


B (o) (t,T)SS 0 ,H 0 ) = 1 


ptl 


£ o K 


F X (X + Y) dX dY 


R < 2 ) (S,i;I 0 ,ti 0 ) = 


E(n) (^? 0 ,-, 0 ) * 


Pi 
i 0 J 
Pi 
i 0 J 


p t 1 


'Ho 


P*) 


F X (X + Y)R (1) (X,Y;| o ,t 1o ) dX dY (2l) 


F X (X + Y)R (n-l) (X,Y;| 0 ,T lo ) dX dY 


(See reference 10.) 
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Except for the fact that F^(0) = co, which causes certain complica- 
tions, the evaluation of the foregoing integrals does not represent any 
essential difficulty. The tabulation of R = R(|,ti;| 0 ,t) 0 ^ is rather 

complicated since R is a function of four variables. 

In appendix A the second term R^ of the expansion (equation (20)) 
is computed and tabulated. 

Remark.— In reference 1, formula (119)^, the expression for a funda- 
mental solution of equation (121) has been derived. The coefficient a 
of the logarithmic term is given by equation (123). It has exactly the 
same structure as the expression for R(| ,11 ; except for the fact 

that the arguments of a are complex quantities, whereas 

are real. In reference 11, Kraft has elaborated a working procedure for 
the evaluation of a which can be used for actual computation of R. 

It becomes in this case simpler since the arguments are real. 


Integral Operator of Type Given by Formula (ll) 

In analogy to the subsonic caBe, where formula (ll) yields a repre- 
sentation for solutions of equation (10) in terms of an arbitrary analytic 
function of one complex variable, there is derived in this section a 
representation distinct from the one" explained in the preceding section 
for solutions of equation (9) (the equation for the stream function in 
the supersonic case) in terms of two arbitrary differentiable functions 
of a real variable. 

Obviously, instead of equation (9), equation (15) for the reduced 
stream function \|r* may be considered. The desired representation can 
be obtained as a consequence of the following theorem. 

Theorem^ .— Let E^ n ^(T) denote functions defined by the following 
recurrence relations: 


E 


< 1 ’<T) --I 


Ht 


f i( t i) dT i 


E 


( n+l) ( T ) = - 


It 

• [ E<n) T lTl ( T l) + S l( T l) B(n) ( T l)] 4T 1 
US 

. . t fS*-., 


> (22) 


Note that in equation (119) and line 10 of page 46 of reference 1, 
by a misprint, a capital .A instead of a, lower case, is used. 

®This theorem has been indicated in reference 3; see chapters 4 
and 5 of this reference. See also reference 12. 
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Here F-^(r) is a given function, the derivatives of which satisfy the 
inequalities 


d n Fx(T) 

dr n 


(n + 1)‘M 


( T o ~ 6 " T ) 


n+2 1 


n = 0,1,2,3, . . . 


(23) 


where e(>0), M = M(e), and t q are conveniently chosen co ns tants , 
Then the series 


00 

= f (l) + E^ n ^(| + 7]) f ( n ) ( £ ) (24) 


n=l 


and 


V ^(t,tl) = g( T| ) + ^ E( n ^(| + T|) g( n ^ ( T]) 

n=l 


( 25 ) 


where f and g are two arbitrary differentiable functions of a real 
variable and 


f [n + l](i) . (Wjy 46i 

Jo 


g 


[ n+ l] (t]) = 


nr\ 

Jo 


g M (^l) dri! 


(26) 


will converge in the intersection of the domains 


-1 < 


T o - ^ - n 


< 1 


and 


-i < — 3 — < i 
T o - 6 - 3 


(27) 


and will represent there solutions of the equation 


H(i|f*) s T|r* + F 1 (g + q)t* = 0 


(28) 
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The proof of this theorem is given in appendix B. As has "been shown in 
reference 13 (see equation (V?)), F^ is given by equation (l6), 

where M, B, and 2A are connected by relations (5) and (17). Thus 
the question of the representation of solutions of equation (15) is 
reduced to the investigation of the behavior of the function F,(t), in 

particular to determining whether its derivatives satisfy the inequali- 
ties (23). This investigation has been carried out in appendix C. If 
in the A, 0— plane the origin is shifted, thus introducing instead of 2A 
the new variable 



and the variables are denoted by and q' as follows. 


I* = (a* + 0) 
V = (a* - 0) 


(29) 


(30) 


then equation (15) is transformed into an equation of the same form 

*Vv + F 2 ^* +' v)+* = 0 (31) 

where F 2 (t), t = |* + q’, satisfies the inequalities (23), for 

F 2 |^2A + — l)l = F^(2A) is an analytic function of A which is 

regular in the strip 

-00 < A. < 00 , 0 < A < jt(h — 1 — l) (32) 


and therefore in particular in the circle of radius t q — e = 

-i - jt — 1^ — e, e > 0. (See fig. 2.) If M = M(e) denotes the 
maximum of F^ ( 2A) in this circle, then by classical results of the 
theory of functions of a complex variable, F^ satisfies the inequali- 
ties (23), where 2A' is substituted for t. Since e can be chosen 
arbitrarily small, the series (equations (24) and (25)) formed for equa- 
tion (31) converge in the domain indicated in figure 3. 

On the other hand, in many instances it is necessary to have repre- 
sentations (equations (24) and (25)) of the solutions which are valid in 
the whole strip. 


— 00 < e < 


0 < 2A < 2jr(i - l) 


(33) 
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One possibility of obtaining such formulas consists in approxi- 
mating F 2 (2A’ ) by a sequence of functions Fg^fSA* ), each of which 

is an entire function of A and therefore (by a classical result of the 
theory of analytic functions of a complex variable) satisfies the 
inequalities (23) and which is chosen in such a manner that 


lim F 2 ^(2A') = F 2 (2A*) 


( 3*0 


in the interval 

0 < 2A < 20 - l)it (35) 

(Concerning another method of representing k = 1,2, see (27), (25) 

in the domain (33), and reference 2, section 5.) Let T|r*( m ) (| » ,tj* ) denote 
the solutions of the equation 


+ F^tt' + tj» )**(“) =0, m = 1,2, 


(36) 


where 


-t (0) < S' < s (0) 


such that ♦*(m)(|,,0) = +*(§•, 0) and ( 0 , r, * ) = ^*(0,^’ ) , 

— t/®) < T) f < As is shown in appendix E, 

lim ♦* (m) (| , ,ri') = **(!’, V) ( 37 ) 

m — > 00 


in the interval 


t (o) (0) (0) , (oy 

-I < I' < i , -11 < V < T1 


(38) 


A sequence of functions F 2 ^ m ^(2A' ) can be obtained in the following 
manner . Let 


p 

F 2 (2A* ) = a Q + a^A* + a^A 1 + . . . 


(39) 
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be the series development of F 2 around the point 2A* = 0 
/ v — 

fi.e., 2A= ^ — ■^ le s® r i e 8 a^A ' 11 obviously converges in 

n =0 

the circle indicated in figure 2 . The 'series 


2 < m) (2A') = ]T 


°n 


n =0 


r(l + n/m) 


A* 


n 


(40) 


will be, for every finite positive number m, an entire function of A* 
and therefore the series (equations (24) and ( 25 )) obtained for equa- 
tion ( 36 ) with F 2 ( m ) given by equation (40) will converge, for 

any m < <», in the strip ( 33 ). On the other hand, according to a 
classical theorem of the theory of functions (see reference 14), the 
relation (equation (34)) holds in this case, so that for the functions 

^♦(m) (| » ,rj» ) obtained in the foregoing manner the relation (equation ( 37 )) 
holds in the whole interval ( 38 ). 

Naturally it is^also possible to approximate F 2 (2a') by using 

different procedures. Some of them are convenient for some special 
purposes, in particular if a problem requires the representation of 
solutions in certain subdomains of the strip ( 33 ). 

SUMMARY REMARKS 


The present paper employs the holograph method for generating flow 
patterns of supersonic flows. The essential feature of this method is 
that the equations are linear and therefore the principle of superposition 
of solutions holds. If, therefore, \|r v (A,6>), V = 1,2,3, . . ., repreV 

sents a set of particular solutions of equation ( 7 ), and A v , arbitrary 
constants, any linear combination 

N 

y A v ^ v (A,e) 

V=1 

is also a solution of equation ( 7 ). E y varying the constants A v , 

v = 1 , 2 , . . ., flows around different shapes can be obtained .1 On the 
other hand, it is often necessary to determine constants Ay to yield 
a flow which approximates that about a prescribed boundary curve, the 
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equation of which is, say, F (x,y) =0. In this case, the constants A v 
can he determined in such a manner that 


M(A) e- J |F[x(A,e;A v ),y(A,0;A v )]| 


d9 


under the condition 


i|r(A,0) = 0 

Here x = x(A^0;A v ) and y = y(A,0;A v ) represent functions of a and 9 
defined in the INTRODUCTION and corresponding to the stream function 
\|r = } 9) in the A, ©—plane. 

As is well known, a solution of the problem does not always exist. 
Under the assumption that \[r v (A,0), v = 1,2, . . ., represents a Bystem 
which is in a certain sense complete, M(A) represents a measure for 
closeness of approximation. In this case the fact that M(A) approaches 
zero when the number N of functions ty v (x,y) increases can be con- 
sidered as a condition for the existence of a solution. 

On the other hand, in many instances, only solutions of the problem 
exist which possess shock lines, and therefore not only the usual condi- 
tion that one of the streamlines coincides with the given boundary must 
be considered but also the condition that these shock conditions are 
satisfied along unknown characteristics. 


Brown University 

; Providence, R.I., November 18, 19^6 
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APPENDIX A 

EVALUATION AND TABULATION OF THE FUNCTION 


The second term R^( in the series for the Eiemann 

function is evaluated in this section. 

F-^ is a function of the variable A = - + ^ and therefore 


F x (x + Ti) dx = r (1) u + T)) - r (1) (£ o + ti) 


P u 

dY 


PI 


pn 


F(X + Y) dX = 


!«■ 


[r (1) (! + y) -r (1) (i Q + Y) 


dY 


l^c 


= r (2) U + ti) - r (2) (? + Ty - r ( 2 ) (i o + ti) 


+ f( 2 Ve~ + 


(*o + %) 


(41) 


where 

1t; l 

F x (t) dr = r^^t-L) -r( 1 ) (t 0 ) (42) 

Jt D 

and 


"Iti 

r ( 1 ) (r) dr = r ( 2 ) (t x ) - r ( 2 ) (t 0 ) (43) 

Jt 
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The fact that F(0) is infinite causes a certain amount of inconvenience 

in the tabulation of the functions and r( 2 t According to equa- 

tion (l 6 ), for k = 1.4, 

F, (B) = ~ °-51 - 0.12B 2 (44) 

1 BP IT B? 

v 

where 

A = 2 \]~6 arc tan — 2 arc tan B (45) 

% 

and. B 2 is given by equation (17), M being the Mach number. 

In table I the values of M, B, and F^ have been tabulated; in 
figure 4, F-|_ is given as a function of M, B, and 2 A. For 
2A = O .3686 the function F-j_ vanishes, that is, Fj [b(0.3686)J = 0; 

this point has been chosen as the initial point for the integration, 
since, as has been pointed out earlier, for the most natural choice 
2A = 0, Fj becomes infinite. 

The values of r^ 1 )(2A) and r( 2 )(2A) for (O. 3686 ) = 

(O.3686) = 0 are given in table II. The interval of tabulation 

of 2A is one four— hundredth. In order to obtain F^ 2 ^ (2A) for inter— 
mediate values ^ the interpolation, formula 

r( 2 )(2A) = r (2) (2A D ) + pD-L + p 2 D 2 (46) 

can be used, where 

r^ 2 )(2A) = r( 2 )(2A Q + ph) = r( 2 )(2A Q ) + ph + |p (p - l)h 2 

= r^)(2A 0 ) + pD-L + p 2 D 2 (47) 

The quantity h is the interval of tabulation at 2 /^ 
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APPENDIX B 


PROOF OF THEOREM GIVEN BY EQUATIONS (22) TO (28) 


The proof is carried out at first under the assumption t ha t the 
series^ (equations ( 2 k ) and (25)) as well as those for their first 
derivatives converge uniformly. (The validity of this assumption is 
proved in this appendix.) A formal computation yields 



n=l 


that is. 


+ ]T f[n]j^.(n) |lT)t 

n=l 

+ ^ n+1 \, + F x E^ n )J (50) 


where is a function of one variable t = 2A f . 

Assuming that the functions are also functions of one 

variable t 9 


(51) 


which yields the foregoing relations, if the assumption is made that 

E( n )(a) = 0 

^In this appendix A,|,t) are replaced by A , ,t , ,T) , . 



r(D = 


= f E 


( 52 ) 
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It remains therefore to prove that the series and "both its formal 
(i.e., term-by-term) derivatives converge uniformly. 

Notation.— If in the interval I there hold for the functions A(t ) 
and 7i(r) and all their derivatives the inequalities: 


|a(t)| <^(t) and 


A(r) will he denoted as a dominant of A(t), which fact is symbolized 
hy writing A « A or 


d n A(T ) 


dr 


n 


<d^(T) 


dr' 


n 


, T€l 


(53) 


A » A, T € I 


(54) 


If e( 1 )(t) 


is given hy 


eW (t) 


n r 


*'l( T l) dT l 
UO 


(55) 


where 


F x (t) » F^t) 


(56) 


then 


II 


IT 

f i( t i) iT i 

n* 

< 


l 

a 

Ja 


*l( T !) dr lldlfr), 


Tel 


(57) 


and also 


Thus 


d n E ( 1 )( T ) 

dT n 


< t€I 
dT n 


(58) 


e( 1 ) (t) « iK 1 ) (t), Tel 


( 59 ) 
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Suppose now that 
nr 


E<n+1)( T ) = 


® (n) Ti T i( Tl ) + fl ( T l) 1<n) ( T i)] dTj + (60) 


where 


E 


( n )(-r) « (t ) f 0 «H( n )(T)^ t € I 


(61) 


Then it follows immediately that 

| E (n+l)( T )| < g(n+l)( T ) 

dE (n+l)( x ) I aE( n+ l)(-r) 


dT 


dT 


(62) 


and, hy considering the corresponding derivatives of [S^ n ^ TT + F^E^ n ^j 
in comparison with [e^ u ^ tt + F-jE^J, it follows that 


E 


(n+l)( T ) « E^ n+ l) ( t) , Tel 


(63) 


which completes the proof "by induction. 
Now hy expression (23), 

F x (t) « F x (t) = 


M 


If H 1 


(n)( T ) 


is given hy 


( T o“ 5 - T ) 


C„Mt 2 


2n 


, T6I 


(6k) 


#n) = n ° 


( T o “ e ~ *) n 


(65) 
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where the C^'s are some conveniently chosen positive constants (to he 

determined later), then an explicit expression may be obtained for E(t ) > 
namely. 


S(n)(T) = CnMT ° g (66) 

( T c - 5 - T ) 

In order to express C n+ j in terms of C n , the right-hand side of 
equation (66) is substituted into equation ( 6 0 ) to give 


E ) ( t ) = C n MT Q ^ 


where 


IT 

n(n+ 1) + C n M+ 0 2 


, _ x n+2 

Ja 

L( T o - e - T ) J 


dT + 


c n + l MT o 2 




.2 


( T o~ €_a ) n C T ° e_T ) 


- * - x \n+2 

(67) 


C n+1 - C n' n + 


MV 

n + 1 , 


( 68 ) 


Thus for sufficiently large values of n 

C n + 1 = C n( n + D (69) 

and therefore for every value of M 

c n +i = M*(n + l) I (70) 

where M* is a conveniently chosen constant. The function f ( £ ) is 
assumed to be differentiable and therefore there exists a constant, 
say, M 2 , such that 


|f(l’)| < Mg for 16*1 g l 2 


Then 


( 71 ) 






in 


= n! 


for U'l = 


( 72 ) 
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Consequently, from equation (24) it is seen that 


V 1 U , ,ti’) « 2M»M2A 0 2 M 


1 + 


U» I g! , 

( T o~ e ~ T ) 


lrl n (n + 1); t 
nl ( T o -e-T ) n 


which converges if 


-1 < 



A» + 6 

( T o - e - 2A») 


< 1 


(73) 


(7*0 


Since e can he chosen arbitrarily small, the series (equation (2k)) 
converges in the domain 


-1 


< A» + 0 
T 0 - 2A* 


< 1 


(75) 


(see fig. 3), where the domain (75) is hounded hy solid lines. 
Similarly, the series (equation (25)) converges in the domain 


-1 


< 


A 1 - e 

t q - 2A 


< 1 


(76) 


which is hounded hy a dashed line so that both series (equations (24) 
and (25 )) converge in the shaded domain. 
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APPENDIX C 

/ 

INVESTIGATION OF THE FUNCTION F 1 (2A) 


This section is devoted to the investigation of the behavior of 
the function F^(2A) in order to determine where its derivatives satisfy 

the inequalities (6k). 


If M varies from 1 to °°, 

and 


B = (M 2 - l) 1 / 2 


varies from 0 to ® 


A = i arc tan (hB) — arc tan B 


(5) 


from 0 to — 1^. Therefore it is necessary to investigate the 
behavior of Fj(2A) on the interval 



(77) 


Since F-^ (given by equation (l6)) is a rational function of B, it is 
sufficient to investigate the function B(2X). 


In order to carry out this investigation it is convenient to continue 
the functions to complex values, so that it will be possible to use 
methods of the theory of analytic functions of a complex variable. 

Instead of equation (5), the complex function 

UP) = X + 1A = | (log 1 ~ £ + i log -- - * ^ (78) ■ 

2 V i + n h h- 1 -?/ 

Z' r»J 

of the complex variable 3 = T + iB will be considered Ip = T, 0 < T < 1 

When P is real it is seen that, if the principal values of the 

logarithmic terms are chosen, then 2 is also real and the foregoing 
equation reduces to: 
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, 1 1 - T 1 .. 

x - 5 l loe m + 5 108 v -i 


h -1 + t' 


A = 0 


h _J - - T 


(79) 


in agreement with the earlier definition of X. On the other hand, if 

0 = iB, B > 0 


B is purely imaginary; then, if the principal values of the logarithmic 
terms are taken once again, there results: 

X = 0, iA = i£— arc tan B + ^ arc tan (hB)J (80) 

which is also in agreement with the previous definition. The question 
which has to he investigated is to determine the domain, of regularity 
of 0(2), inverse to equation (77). 

A classical theorem concerning the inversion of ana lytic functions 
states that, if w(z) is a function of the complex variable z which 
is regular at z = z Q and has a nonzero derivative there, then, in some 

sufficiently small neighborhood of the point w = w Q . = w(z Q ), it is 

possible to invert the function w(z), that is, to express z as a 
regular function of w. (See reference 15 , p.^lk2. ) At every point 
except 0 = ±1 and 0 = ±h, the function 2(0) is regular, so that 
th 9 zero points of the derivative d2/d0 have to be determined. 

A formal computation yields 


d2 

dp 


1 

2 


1 
+ — 


1-? i + 0 Hh-i + 




0 h~ 1 - p 


(l - h 2 )^ 

(l - 0 2 )(l - h^p 2 ) 


( 81 ) 


Obviously, 


tc = 0 for f-f = 0 and 0 = « ( 82 ) 

d0 

The discussion would become quite complicated if all^the possible 
branches of the function 2(0) and of the inverse 0(2) were to be 
considered, because each summand of equation (8l) is an infinitely 
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many— valued function of 3, having logarithmic singularities at ±1 

and +h -1 , respectively. It is therefore convenient to render 7(3) 
single-valued by the artifice of starting at 3=0 with the principal 

values of log(l — 3), log(l +3), log(h — 1 +3) and log(h — 1 — 3) 
(i.e., with the real values of the logarithms) and slitting the complex 
3— plane from 1 to +°° and from —1 to — along the real axis . 

Since h — ^ > 1, it is clear that the slit plane contains no branch 
points and therefore 7 is uniquely defined by beginning with the 
aforementioned determination of 7(0). 

A/ 

Now, a direct investigation shows that the 3— plane, slit in the 
manner Just described, is mapped into the band 

— 00 < X < «>, — (h~ 1 — l)jt < A < (h — 1 — l)n (83) 

However, the mapping is not bi-unique, for at the branch point 3=0 
the mapping is easily seen to be 3 to 1. (See appendix D.) Thus, the 
simply covered slit 3— plane is mapped into the simply covered band Just 
defined. The only singularities of ”3(7) in the strip 


— 00 < X < 00, 0 < A< (h — l) jt (8h) 

are the^jpoints 7 = 0 where 3(7) has a branch point end 7 = (h — — l)ni 
where 3(7) has a pole. 
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APPENDIX D 


THE DEVELOPMENT OF 0(2) AT THE ORIGIN 2 = 0 


The simplest method of proving that 2 = 0 is a "branch point of $(2) 
of third order is to determine the series development of this function 
at 2=0. In order to obtain this development, the principal branch 
of 2(0), namely. 




arc tan (0/i) arc tan (0h/i)l ( 85 ) 


has to be developed in a power series to give 

2 = - £ (l - h 2 ) - 2 (1 _ h 1 *) - 21 (l - h 6 ) . 
3 5 7 


1 - h c 


.33 3 (l-hM p 3 (l-h 6 ) ?7 

5 (l - h 2 ) ^ Ml-h 2 ) 


Introducing a new variable s. 


s^ 2, that is, s = ( — 

1 - h 2 V 1 - 


yields the development 


8 = C ^0 + C ^- 3 + C^fP + . 


Since c as is easily seen, is not zero, it follows from the 
general theory of analytic functions that this series can be inverted: 


$ = djS + d 2 s^ + d^s^ + . . 
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In particular, for h = l/\| 6 , 

C 1 = c 3 = = 1, d^ = — 0.2333 • • (9l) 


\ 


10 Th 3 values of c n and d n , n = 1 , 3 , . . 
puted "by the author and will he published elsewhere. 


25 have been com- 


NACA TN No. 1875 


37 


APPENDIX E 

APPROXIMATIONS OF SOLUTIONS 


In the present report , solutions of equation (15) have heen approxi- 
mated at several places by solutions of an equation in which the coeffi- 
cient F^ is replaced by another coefficient which approximates it. The 

legitimacy of such a procedure is established by the following theorem. 
Theorem .— Let i|X^(i,T|), k = 1,2, be solutions of the equations 

+ f( k )i|r( k > = 0 . (92) 


where the f^'s are given functions of | and t] which satisfy 
throughout the rectangle 0 < | < a, 0 <_ q < b the inequalities 


| f (l) _ f (2) | < e ^ |f(k)|< m (93) 

( V ) 

€ and m being positive constants. The functions i|r v • are also to 
satisfy the conditions: 

i (l) (!,0) = * (2) (|,0) = X x (5) for 0<l<a (9V) 

and ' 

(0,q) = ^(ti) for 0 < TJ < b (95) 


where X^ 

intervals 
tion that 



and Xg sure continuous functions prescribed on the closed 

(0,a) and (0,b), respectively, subject only to the condi— 
X^O) =X 2 (0). Let x a be a positive number such that 
| Xg j <<x. Then there exists a positive number M, depending 


only on e, a, m, and ab (the area of the given rectangle) such t ha t 
for fixed values of a, m, and ab, M approaches a finite (positive) 
limit as € approaches zero, and such that, throughout the given rectangle 
the following inequality holds: 



IRroof According to the general theory of linear hyperbolic partial differential equations, it 

is known that the functions satisfying expressions (92), ( 93 ), and ( 95 ) ina y he expressed in 

the form of the following infinite series, which is valid throughout the given rectangle: 
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so that, according to expression (93), |8(S,t))| < e. Now set k equal to 1 and 2, respectively, 

in equation (97) and subtract the two equations so obtained (pairing off corresponding terms and 

replacing f^ 1 ) by f( 2 ) + 6 ). The result obtained is: 




(TQT) T ^P T 5P S ^P S IP ( S U' S s)d(V 5 s) 
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And it is clear how a s imi lar reduction can "be effected in all the 
succeed ing terms, so that no sums appear in any of the integrands. 

Now, if each term in the series on the right-hand side is replaced 
by its absolute value, an upper bound is obtained for the absolute value 
of the lefir-hand side of equation (100); and this inequality continues 

to hold if the functions p, and f^ are replaced by their respec- 

tive upper bounds 3a, e, and m. (The upper bound for p follows 
directly from equation (98) and the definition of a.) Therefore, from 

equation (100) there is obtained the inequality: 

* 

LW l< 3ae {^3a(6V2€m)i 2 3 g + 3a(e 3 - , -3€ g m- | -3€m g )lV +j _ . (102 ) 

( 2!) 2 ( 3!) 2 


It is not difficult to see that the nth term of the right-hand side 
of equation (102) is simply: 


3a [(e + m) n - mJ| n T} n 


(n!) 


2 


( 103 ) 


The inequality (102) holds a fortiori if | and t] are replaced 
by their respective upper bounds a and b. Thus, there is finally 
obtained for all points in the given rectangle, the inequality 


Tl/ 1 ) _ 


n=l 


3a[(c + m) n — m n ] (ab) n 


(n!)‘ 


( 104 ) 


This series is easily shown, by elementary tests, to converge. Further- 
more, it is clear that the quantity e can be factored out of the 

expression [(€ + m) n — nf 1 ], so that the series may be written as the 
product of e by a new series, which is again easily shown to converge 
for all values of e. Designating this new series by M, the conclusion 
of the theorem is obtained. 


In less precise language than that used in the formulation of the 

theorem, the result obtained may be expressed by saying that, if the 
coefficients f (4) of the two differential equations (equation (92)) 
approximate each other very closely in a certain region and if the 
respective solutions satisfy the same side conditions (equations (94) 
and (95) ), then the solutions approximate each other very closely 
throughout the given region. 
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In particular it follows from this theorem that, if the sequence of 
coefficients of the given hyperbolic equation converges uniformly, the 
corresponding sequence of solutions also converges uni formly to the 
solution of the limiting equation. 
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TABLE II.- F-^ H 1 ), AND H 2 ) AS FUNCTIONS OF 2 A, M, AND B 


2 A 

M 

B 

F 1 

r (l) 

r (2) 

0.00 

1.0000 

0 

00 

— 00 

— 00 

.04 

1.0899 

. 433 ^ 

84.078 

- 2.6449 

-0.1211 

.08 

1.1460 

.5597 

19.833 

-.'9974 

-.0561 

.12 

1.1951 

.6544 

8.1146 

-.4859 

-.0279 

.16 

1.2406 

• 734 l’ 

4.0787 

-.2537 

-.0136 

.20 

I.2838 

.8051 

2.2408 

-.1315 

-.0061 

. 2 b 

1.3255 

.8701 

1.2579 

-.0634 

-.0024 

.28 

1.3663 

.9310 

.6737 

-.0257 

-.0007 

.32 

1.4063 

.9887 

.2989 

-.0067 

-.0001 

.36 

1.4457 

1.0441 

.0441 

-.0002 

.0000 

.bo 

1.4848 

. 

1.0976 

-.1375 

-.0022 

-.0000 

M 

1 . 533 ^ 

1.1625 

-.3003 

-.0134 

-.0004 

.50 

1.5819 

1.2257 

-.4182 

-.0315 

-.0015 

.55 

1.6303 

1.2876 

-.5077 

. -.0548 

-.0036 

.60 

1.6789 

1.3486 

-.5784 

-.0820 

-.0070 

.65 

1.7277 

1.4089 

-.6365 

-.1124 

-.0119 

.70 

1.7768 

1.4687 . 

— .6860 

-.1455 

-.0183 

.75 

1.8264 

1.5283 

-.7294 

-.1809 

-.0265 

.80 

1.8765 

1.5878 

-.7686 

-.2184 

-.0364 

.85 

1.9271 

1.6474 

-.8050 

-.2577 

-.0483 

.30 

1.9784 

1.7071 

-.8394 

-.2988 

-.0622 
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I _l 1 1 1 1 : 1 ► B 

0 . 663 . 922 1.14 1.35 1.54 1.73 


Figure 4.- F^ as a function of M, B, and 2a 


